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The complete low-lying positive charge conjugation glueball spectrum is obtained 
from QCD. The formalism relies on the construction of an efficient quasiparticle gluon 
basis for Hamiltonian QCD in Coulomb gauge. The resulting rapidly convergent Fock 
space expansion is exploited to derive quenched low-lying glueball masses with no free 
parameters which are in remarkable agreement with lattice gauge theory. 



Introduction The scalar glueball has been called the fundamental particle of 
QCD[1.. Indeed, its existence and nonzero mass are a direct consequence of the 
nonabelian nature of QCD and the confinement phenomenon. It is clear that finding 
and understanding the scalar (and other) glueballs is a vital step in mastering low 
energy QCD. 

Recently quenched lattice gauge theory has been able to determine the low lying 
glueball spectrum with reasonable accuracy [2| (only very preliminary determinations 
of other matrix elements have been attempted). These data serve as a useful bench- 
mark in the development of a qualitative model of the emergent properties of low 
energy QCD. The models may then be used to guide experimental glueball searches. 

Previous models of glueballs have relied on ad hoc effective degrees of freedom such 
as flux tubes |3] , bags jl] , or constituent gluons |5] . We note that some of the models 
listed in Ref. construct states with massive gluons, while others either use trans- 
verse gluons or dynamically generated gluons masses. Models in the former category 
contain spurious states due to the presence of unphysical longitudinal gluon modes. 
Sum rule computations of glueball properties exist however, they are based on 
phenomenological properties of the spectrum. Finally, the conjectured duality be- 
tween supergravity and large-N gauge theories has been used to compute the glueball 
spectrum in non-supersymmetric Yang-Mills theory by solving the supergravity wave 
equations in a black hole geometry [7j. Unfortunately all of these approaches suffer 
from weak or conjectured connections to QCD. 

We present a computation of the positive charge conjugation glueball spectrum 
which arises from QCD and is systematically improvable. The computation is based 
on the formalism presented in Ref. || in which the QCD Hamiltonian in Coulomb 
gauge is employed as a starting point. Coulomb gauge is efficacious for the study of 
bound states because all degrees of freedom are physical (there are no ghost fields 
in this gauge) and a positive definite norm exists 0. Furthermore, resolving the 
Coulomb gauge constraint produces an instantaneous interaction (the nonabelian 
analogue of the Coulomb interaction) which, as shown in Ref. jH] may be used to 
generate bound states. Because the temporal component of the vector potential is 
renormalization group invariant in Coulomb gauge (this is not true in other gauges), 
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the instantaneous potential does not depend on the ultraviolet regulator or the renor- 
malization scale • This fact permits a physical interpretation of the instantaneous 
potential which is a central aspect of our formalism. 

The pure gauge QCD Hamiltonian may be written as P] Hqcd — H + 5H with 



H = ~ I r/x 



E? 2 + B' 



+ i| C /xrfyp a (x)K(°)(x-y)p a (y) 



and 



5H = V 3g + V ig +Vj + V c (2) 

Here B = V x A is the abelian part of the chromomagnetic field and E = —d/dtA 
is the chromoelectric field. The third term in Hq represents the nonabelian, instan- 
taneous Coulomb interaction between color charges, p a = — f abc ~Ei b ■ A c , mediated 
by an effective potential K° computed by taking a vacuum expectation value of the 
Coulomb kernel, 

iT°(x - y)S ab = g 2 {*\ [(V • J D)- 1 (-V 2 )(V • D)- 1 } |*>. (3) 

L J x,a;y,o 

with D ab = <5 afe V — gf abc A c being the covariant derivative in the adjoint representa- 
tion. For the vacuum wave functional, ^[A] = (A\^f) we take a variational ansatz, 

1 f d 3 k 

with the variational parameter u(k) determined by minimizing the vev of H. The 
correction terms in 6H include V 3g and V± g which are the three and four gluon opera- 
tors originating from the difference between the full and the abelian chromomagnetic 
field. Vj denotes a contribution from the Faddeev-Popov determinant in the kinetic 
term. The effects of Vj and the Faddeev-Popov determinant in the functional inte- 
grals have recently being studied in Ref. ^U] where it was found that its effects can 
be effectively included in the variational parameter u(k). Finally, Vq is the difference 
between the Coulomb operator and its vev, K°. In the calculation of the glueball 
spectrum it results in operators mixing two and three, quasi-particle wave functions. 
We note that after renormalization the coupling g is absorbed into the Faddeev-Popov 
operator, which then defines the Coulomb gauge analog of a ghost propagator [HlllUj. 



9[A] = exp J -|^A a (kV(fc)A a (-k)) (4) 
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The renormalized effective potential K° is fixed by comparing with the quenched lat- 
tice QCD static potential. A very accurate representation of the static confinement 
potential is achieved jS]. 

The variational vacuum defined above also specifies a Fock space of quasiparti- 
cle excitations corresponding to effective gluonic degrees of freedom (which we call 
quasigluons) . Such quasigluons obey "massive" dispersion relation in the variational 
vacuum and therefore improve the description of gluonic bound states since mix- 
ing between states with different number of quasiparticles is suppressed due to their 
effective mass. 

The calculation of the vev of the Hamiltonian and the properties of the quasi- 
particle excitations were discussed in Refs. [El HH EI] ■ These require solving a set of 
coupled integral Dyson equations and as a result one finds that the function u(k), 
which in a free theory is given by u(k) = k, becomes finite as k — > 0. The value o>(0) 
can be related to the slope of the static potential at large distances. 

Fock Space Expansion and the Glueball Spectrum. The quasigluons which emerge 
in the analysis of Ref. [S] set the QCD scale parameter via the low momentum disper- 
sion relation r u(k — > 0) = 1.4, where r is the lattice Sommer parameter. Using the 
Regge string tension or p mass to fix the scale then gives uj(0) = 600 — 650 MeV. It 
is natural to interpret this scale as a dynamical gluon mass ^3]. Thus the formalism 
of Ref. jH] provides a justification of a Fock space expansion in terms of quasigluons 
and gives the leading instantaneous interaction between the quasigluons. In view of 
this it is natural to attempt a description of low lying glueballs in the pure gauge 
sector of QCD. 

In this approach positive charge conjugation glueballs are dominated by the two 
quasigluon contribution. These may mix with three and higher quasigluon states via 
transverse gluon exchange (and, in general, via any term in SH). Mixing with single 
quasigluon states is excluded because color nonsinglet states are removed from the 
spectrum due to infrared divergences in the color non-singlet spectrum [TT]. Finally, 
the scalar glueball is orthogonal to the vacuum due to the form of the gap equation. 

The resulting bound state equations are shown in Eq. (jSJ). There is one orbital 
component of glueball wave function for J p = + and two for J p = (even > 2) + . 
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These are denoted by ipi(k), i — 1,2. The first term on the right hand side of Eq. (J5J) 
represents the quasigluon kinetic energy (the gluon gap equation has been employed 
to simplify the expression), the second term is the quasigluon self energy, and third 
represents the interaction between quasigluons in the channel of interest. 



y 



with 



r k 2 dk , ., . „ N c ^ f 

N c 



k 2 dk q 2 dq uj(k) 



k 2 dk q 2 dq (u(k) + u(q)) 2 
(2tt) 3 (2tt) 3 u{k)uj(q) 



(2tt) 3 (2tt) 3 uj(q) 
^(q)K tl (q,k)^(k) 



4 2 
- 3 Vo + - 3 V, 



E 



k 2 dk 
(X) 3 



3J 2 + 3 J - 2 
(2J- l)(2J + 3) 

(J+l)(J + 2) 
2(2J + 3)(2J+1 



Vj + 



J(J-l) 



2(2J- 1)(2J+ 1) 



V.. 



J-2 



(5) 



(6) 



and 



K 
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+ 



3(J + 2)(J-1) 
(2J-l)(2J + 3) 

2(2J+ l)(2J + 3) 



Vj 



(J + 2)(J+1) 
2(2J + 1)(2J- 1 



-V, 



J-2 



V 



J+2 



(7) 



K 



12 



^21 = v/(J-l)J(J + l)(J + 2) 

1 Vj+2 + 



2(2J + 3)(2J+ 1 
1 



2(2J+ 1)(2J-1 



-Vj- 



-V, 



(2J + 3)(2J- 1) 

The bound state equations for other glueballs are as in Eq. with the exception 
that the wave function index takes on a single value. For these cases the interaction 
kernels are as follows: 

J p = (odd > 3) + (there is no 1 + gg glueball) 



K 



J+2 
2J+ 1 



V, 



j-i 



J - 
27 



-v 



J+l 



(9) 



5 



J p 



{even > 0) 



K 



J 



J+i 

2J+ 1 



(10) 



2J+ 1 




(11) 



In all these relations the interaction is defined as 



V L (q,k) = 27T [ dxK ( - 0) (q,k,x = k-^)P L (x) 



(12) 



and the potential K^°\q,k,x) is that derived in Ref. [S] with the QCD scale chosen 
to be uj(0) = 600 MeV. Finally we note that there are no J p = (odd)~ , or C = — 
glueballs at lowest order in the Fock space expansion. 

Higher Order Corrections. It is of course desirable to test the efficacy of the Fock 
space expansion employed here by explicitly checking the effect of coupling to the three 
or higher quasigluon spectrum. This is a difficult coupled channel problem and we 
therefore content ourselves with a perturbative evaluation of these effects in this initial 
study. Specifically, the energy shift 5E n = £ m \{gg\SH\m{ggg))\ 2 / \E n - E m ) must 
be evaluated. Duality implies that when the energy transfer is large, (E n — E rn ) > 
A, where A is of the order of the QCD scale, this sum may be evaluated in its 
perturbative form (with partonic gluons in the intermediate state). We compute here 
the effects of the three-gluon coupling from SH. This is the leading interaction in 
terms of expansion in the coupling constant. After renormalization g 2 /4:7c — > a(p 2 ) 
where p 2 represents a characteristic momentum in integrals when computing matrix 
elements. The running coupling expansion for the remaining sum over low energy 
modes is certainly less justifiable, however as shown in examples in Ref. jS], such soft 
corrections also seem to be small. For numerical efficiency the low energy part of the 
guasigluon exchange was approximated with a local four-gluon interaction (we note 
that this effective interaction also accounts for the four-gluon interaction present in 
the Hamiltonian) 




d 3 k\ d 3 k 2 d 3 k 3 d 3 k4 
(2tt) 3 (2tt) 3 (2tt) 3 (2tt) j 



^(27r) 3 5(k 1 + k 2 + k3 + k 4 )- 
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exp (-(*? + k 2 2 + k\ + fc 4 2 )/A 2 ) ^{(kO^Cka)^ (k 3 )^ e ( k 4)- 

(13) 

where C is a dimensionless parameter of the order of g 2 (A) ~ 1. Standard effective 
field theory techniques were subsequently employed: the factorization scale A was 
chosen and the coupling C was fixed by comparison to the lattice scalar glueball 
mass. Other mass predictions then follow. The scale A was then varied to ensure 
that the procedure is stable and that the coupling remains 'natural' (of order unity). 
To maintain consistency the effect of these terms should be incorporated into the 
quasigluon gap equation and the gluon self energy. However, it may be shown that 
the effect of contact terms are canceled in the bound state equation when the gap 
equation is used to simplify the quasigluon kinetic and self energies. This is not true 
for the high-momentum gluon exchange terms which add a UV dominated correction 
to the single gluon kinetic energy. These have negligible effect on low energy spectrum 
and have been ignored. 

Results and Conclusions. The lowest order predictions for the quenched positive 
charge conjugation glueball spectrum are presented in Table 1. We stress that there 
are no free parameters in this computation; the renormalization group parameters and 
the scale were fixed by comparison to the Wilson loop static interaction [S]. Although 
one may anticipate splittings on the order of 100 MeV due to coupled channel effects, 
the general agreement with lattice data is quite good (the x 2 P er degree of freedom 
for the six measured lowest spin-parity states is 1.5). Nevertheless we note that the 
authors of Ref. j2] state that the 3 ++ may have significant mixing with higher states 
and the quoted 4 ++ mass should be regarded as preliminary. 

Although it appears to be difficult to push lattice mass computations above 4 GeV 
it would be interesting to measure the quenched 4~ + glueball mass to test the pre- 
diction made in Table 1. Lastly we note that all radial excitations lie roughly 1 GeV 
above their respective ground states, except the 4 ++ . We have no explanation for 
this curiosity, but note that it implies that lattice extractions of the 4 ++ mass must 
be made with great care. 
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Table 1: Glueball Spectrum 



state 


this work 


LGT 


Ref 




inn mi "Yin en 


(GeV) 






0++ 


1.98 


1.73(5)(8) 





1 






1 754(851(86) 


17J 






1 627(83) 


m 






1 686(24V10) 


2 


11 






1.645(50) 


2 


a 






1 61(71(13) 


1 


3 


u 


3 9£ 




2 




2++ 


2.42 


2 40(2 51(12) 

9 ZL1 7l^fiVl 1 7^ 

2.26(12)(18) 

2.380(67)(14) 

2.337(100) 


[2] 

en 

CHI 
[E3 

[2U 

H3 


2+ +/ 


3.11 


3.499(43) (35) 


[E! 


o-+ 


2.22 


2.59(4)(13) 
2.19(26)(18) 



CHI 


o- +/ 


3.43 


3.64(6)(18) 


2 




2-+ 


3.09 


3.10(3)(15) 


V 




2- +/ 


4.13 


3.89(4)(19) 
3.69(4)(18) 2 


V 




3++ 


3.33 


V 




3++/ 


4.29 








4 ++ 


3.99 


3.65(6)(18) 


HQ 


4++' 


4.28 








4 -+ 


4.27 








4-+' 


4.98 
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We note that the degeneracy between parity states reported in the first of Refs. [3] 
is not seen here. We suspect that the degeneracy is due to the nonrelativistic expan- 
sion of the interaction kernel made in that reference. 

As stated above, coupled channel effects are expected to modify the spectrum at 
the 10 — 100 MeV level. As an initial estimate of the size of these effects we simply set 
C = (Eq. IT3|) and evaluate the energy shift due to perturbative one gluon exchange. 
As expected, the scalar glueball mass experiences the largest shift, with a reduction in 
mass of roughly 200 MeV (from 1980 to 1790 MeV). It is gratifying that this brings 
the scalar glueball into excellent agreement with lattice gauge theory. We proceed 
by incorporating the effective contact interaction. The factorization scale was varied 
between 0.25 and 10 GeV, stable results were found between 1 and 8 GeV, with a 
value of C given roughly by —0.4 in this range. We find that the tensor glueball mass is 
reduced by roughly lOOMeV, while other masses experience somewhat smaller shifts. 
Thus it appears that low lying glueballs are indeed dominated by their two quasigluon 
Fock components. However it is clear that a careful examination of coupled channel 
effects and better lattice data are required to make a definite statement about the 
efficacy of our approach. 

Fluctuations of the topological charge density have pseudoscalar quantum num- 
bers. This raises the possibility that the QCD anomaly affects the lightest h glue- 
ball mass. Topological effects have so far not been incorporated into our formalism. 
Doing so would require modification of the vacuum ansatz to reflect the identification 
of gauge equivalent field configurations at the boundary of the fundamental modular 
region. This allows contributions from field configurations with nonzero topologi- 
cal charge. Indeed a cross-over between the _+ and 2 ++ glueball masses has been 
observed on the lattice as a function of the renormalised coupling [22] if boundary 
conditions are imposed. 

Further aspects of the gluonic structure laid out in Ref. jH] may be investigated 
by an examination of the adiabatic potential surfaces of heavy quark hybrids (this 
probes nonperturbative gluon-confinement potential couplings). Extensions to the 
light hybrid spectrum will prove of interest to searches for these new states at Jef- 
ferson Lab. Finally, the short range structure of the meson sector is dominated by 
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coupled channel effects and nonperturbative gluodynamics. The wealth of experi- 
mental information in this sector will provide a definitive test of the dynamics being 
advocated in our approach. 
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